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In the recent paper we explained why the maximum bulk resistivity of topological insulators (Tls)is 
so smalffl Using the model of completely compensated semiconductor we showed that when Fermi 
level is pinned in the middle of the gap the activation energy of resistivity A = 0.3(E g /2), where E g 
is the semiconductor gap. In this paper, we consider strongly compensated n-type semiconductor. 
We find position of the Fermi level p calculated from the bottom of the conduction band and 
the activation energy of the resistivity A as a function of compensation K, and show that A = 
0.3(E C — p) holds at any < 1. At the same time Peltier energy (heat) II is even smaller: 

n ~ 0.5 A = 0.15(E C — p). We also show that at low temperatures the activated conductivity crosses 
over to variable range hopping (VRH) and find the characteristic temperature of VRH, Tes, as a 
function of 1 — K. 



I. INTRODUCTION 

The three-dimensional (3D) topological insulator 
(TI^PHU has gapless surface states, which host a spectrum 
of quantum transport phenomena^! While a number of 
crystals have been identified to be 3D TIs, most of them 
are poor insulators and the bulk of TI crystals of sub- 
stantial size (> 10/im) shunts the surface conductivity. 
The current literature^"^ broadly discusses how one can 
achieve a bulk-insulating state. 

Typically as-grown TI crystals are heavily doped re- 
type semiconductors. To make them insulating they are 
compensated by acceptors. With increasing compensa- 
tion K = Na/Nu, where Nr> and Na are the concen- 
trations of monovalent donors and acceptors, the Fermi 
level shifts from the conduction band to inside the gap 
and then into the valence band. When compensation 
of donors is complete, K = 1, the Fermi level is in the 
middle of the gap and the most insulating state of TI is 
reached. The hope is that for a TI with the gap E g ~ 0.3 
eV the resistivity obeys the acivation law 

p = po cxp(A/fc B T) (1) 

with activation energy A = E g /2 ~ 0.15 eV, so that at 
room temperatures and below TI is well insulating. 

The typical experimental situation at K = 1 is however 
frustrating^. In the range of temperatures between 100 
K and 300 K the resistivity is activated, but with three 
times smaller than expected activation energy A ~ 50 
meV. At T ~ 100 K the activated transport is replaced by 
variable range hopping (VRH) and resistivity grows even 
slower with decreasing T. Finally, at even smlaller T < 
50 K resistivity saturates around p(T) < lOfkm. This 
means that in spite of complete compensation even at 
helium temperatures conductance of TI samples thicker 
than 10pm is dominated by the bulk. 

In the recent papei^ we suggest an explanation of 
anomalously large bulk conductivity of TI at K = 1. 
We assume that both donors and acceptors are shallow 
and randomly positioned in space, and use theory of a 



completely compensated semiconductor (CCS)M 1 9'. This 
theory is based on the idea that at K = 1, when al- 
most all donors and acceptors are charged, random spa- 
tial fluctuations of local concentration of impurities re- 
sult in large fluctuations of charge. Their potential is 
poorly screened, because of vanishing average concentra- 
tion n — Nd — Na of electrons. Huge fluctuations of 
the random potential bend conduction and valence band 
edges and in some places bring them to the Fermi level, 
creating electron and hole puddles, which in turn non- 
linearly screen the random potential. As a result, the 
amplitude of potential fluctuations is limited by E g /2 so 
that the ground state of completely compensated semi- 
conduct or sho wn in Fig. [T] reminds a network of p-n 
junction d 18 l 19 [ Characteristic size of these p-n junctions 
R for E g ~ 0.3 eV, N D = 10 19 cm" 3 , the dielectric 
constant k = 30 is R = 150 nm > N D 1/3 = 4.6 ran 1 , 
i.e. we deal with a very long range potential. As a re- 
sult, the resistivity can be dramatically different from 
the one for the flat bands picture of TIM First, at rel- 
atively high temperatures activated conductivity is due 
to electrons and holes activated from the Fermi level to 
their corresponding classical percolation levels (classical 
mobility edges), E e and E^, in the conduction and the va- 
lence bands. According to numerical modelling at K = 1 
1 the activation energy A ~ 0.15.E ff , because E e and 
Eh are substantially closer to the Fermi level p than the 
unperturbed by random potential bottom of conduction 
band E c and ceiling of the valence band E v (Fig. [T^, ). 
Second, at low enough temperatures electrons and holes 
can hop (tunnel) between puddles, so that variable range 
hopping conductivity replaces activated transport. We 
showed that the activated resistivity crosses over directly 
to Efros-Shklovskii (ES) law^ 

p = p Q eMT ES /T) 1 '\ (2) 

where Tes = Ce 2 /fc,BK£, e is the electron charge, £ is 
localization length of the states with the Fermi level en- 
ergy and C = 4.4 is a numerical coefficient. Together 
our results for the activated and VRH resistivity estab- 
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lished the universal upper limit of resistivity p(T) one can 
achieve for a 3D TI compensated by shallow inpurities. 





FIG. 1. Energy diagram of a) completely compensated semi- 
conductor (K — 1) and b) stronly compensated semiconduc- 
tor (1 — K <C 1) with the gap E g . The upper and the lower 
straight lines indicate the unperturb positions of the bottom 
of the conduction band, E c , and the ceiling of the valence 
band E v ; the middle straight line corresponds to the Fermi 
level p. Meandering lines represent the band edges, which are 
modulated by the fluctuating potential of charged impurities, 
R is characteristic size of potential fluctuations. Percolation 
levels for electrons, E e and holes, Eh are shown by dashed 
lines. Puddles occupied by carriers are shaded. Shallow im- 
purities levels are not shown because they practically merge 
with band edges. 

In this paper, we change our focus from maximum pos- 
sible bulk resistivity of completely compensated semi- 
conductor at K = 1 to the more practical question of 
dependence of bulk resistivity of strongly compensated 
semiconductor (SCS) on A at 1 — A' <C 1. Indeed, with 
existing methods of growth of TI samples one can not 
get K = 1 exactly and it is important to know how sta- 
ble K = 1 resistivity results are in the case 1 — K <C 1. 
For example, one can ask at which 1 — K activation en- 
ergy A is twice smaller than at K = 1. For definitcncss 
we consider n-type SCS, where concentration of electrons 
n = Nd — Na <C N d and 1 - K < 1. We model numeri- 
cally ground state of such SCS and its resistivity using al- 
gorithms similar to Refffl We find that in agreement with 
analytic theory^ when 1 — K grows screening of random 



potential improves and its correlation length R decreases. 
The amplitude of the random potential decreases as well, 
hole puddles shrink and eventually vanish and the chem- 
ical potential p moves up, so that E c — p decreases. One 
can say that with increasing (1 — A") screening happens 
by bending of the conduction band only while all accep- 
tors remain occupied by electrons and negative. All these 
changes are illustrated by transition from a) to b) in Fig. 

HI 

As a result of these changes with growing 1 — K the 
activation energy A decreases. We find that the relation 
A = 0.3(A C - p) obtained in Ref. 1 for K = 1 remains 
valid for 1 — K <C 1 (see Fig. [7] below) as well. (In p-type 
semiconductor where A — Nd/Na similar relationship 
A = 0.3(/i— E v ) takes place.) By A = 0.97 the activation 
energy A is at least two times smaller than at A = 1. 
This result shows that achieving maximum resitivity with 
A = 0.15E g is problematic. It also explain the origin of 
large scatter of magnitude of A among TI samples^. 

Our prediction that A = 0.3(A C — fi) in principle can 
be directly compared with experiments in TIs. Indeed, 
for each A the position of the Fermi level (E c — p) can 
be found via measurements of the surface concentration 
of electrons in the surface gapless state using Shubnikov- 
de-Haas oscillations. 

At low temperatures we find numerically direct cross- 
over from acivation to ES VRH. We also find how Tes 
being correlated with A decreases with 1 — K. 

Our assumption of random distribution of impuruities 
is crucial for this theory. Usually for samples made by 
cooling from melt distribution of impurities in space is 
a snapshot of the distribution, which impurities have at 
higher temperature, when diffusion of impurities prac- 
tically freezes. In semiconductors with a narrow enough 
gap at this temperature there is a concentration of intrin- 
sic carriers larger than the concentration of impurities. 
Intrinsic carriers screen Coulomb interaction between im- 
purities, so that impurities remain randomly distributed 
in space. At lower temperatures when intrinsic carri ers 
recombine impurities are left in random position d 18 l 21 l If 
diffusion freezes at T ~ 1000A it is reasonable to assume 
that impurities are randomly positioned in a semicon- 
ductor with Eg < 0.3 eV. This justifies the use of this 
theory for typical TIs. Our results are applicable to other 
narrow gap semiconductors, for example, to InSb. His- 
torically, large effort was made to make InSb insulating 
via strong compensation. The goal was to improve char- 
acteristics of InSb based photo-detectors. Results were 
again frustrating: the dark resistivity was too small. Our 
results are in a reasonab le agr eement with transport ex- 
perimental data for InSrJ 22 | 23 [ 

The plan of the paper is as follows. In section II we 
formulate the model, explain the algorithm of mumerical 
simulation of the pseudo-ground state and present results 
for the density of states (DOS). In section III we present 
our algorithm for computation of hopping conductivity, 
analize our results and arrive at a small activation energy 
of conduction band resistivity A = 0.3(A C — fi). We aslo 
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evaluate the localization length of states close to Fermi 
energy and estimate characteristic temperature of ES law 
T E s- In Sec. IV we estimate the thermopower of strongly 
compensated semiconductor and show that Peltier energy 
(heat) n ~ A/2 = 0.15(£ c - /jl). 



II. THE MODEL, PSEUDOGROUND STATES 
AND THE DENSITY OF STATES 

To model a heavily doped SCS we create a cube filled 
by 20000 donors and 20000-Kf acceptors that are ran- 
domly positioned in space. We numerate all donors and 
acceptors by index i and use n$ = or 1 for number of 
electrons residing on a donor or an acceptor. In addi- 
tion we use a variable /, to discriminate between donors 
(fi = 1) and acceptors (/, = —1). The Hamiltonian of 
our system is 

H = E ^ •/'<"< + E vim)** ■ ( 3 ) 

where qi — (/$/2 — n,- + l/2) is the net charge of site i and 
all energies are defined relative to the Fermi level. The 
first term contains the difference between the energies of 
shallow donors and acceptors which for the case of shal- 
low impurities is very close to the semiconductor gap E g . 
The second term of H are the sum of interaction energies 
of charged impurities. If two impurities are at distance 
r » as, where as is the Bohr radius of impurity states, 
one can use the Coulomb interaction V(r) = e 2 /nr. For 
a pair of empty donors, one donor shifts down the energy 
of the electron on the other by energy V(r) — —e 2 /nr. 
This classical form for V(r) is good for a lightly doped 

— 1 /3 

SCS. But in a heavily doped SCS, where as > N D 
most impurities have at least one neighbor at distance 
r < flg and quantum mechanical averaging over electron 
wave function becomes important. (This is why uncom- 
pensated heavily doped semiconductor is a good metal). 
For example, such a donor pair cannot create a state 
deeper than one of the helium-like ion with binding en- 
ergy 4Eb, where Eb — e 2 /2kclb is the binding energy of 
the shallow donor state. Here we deal with heavily doped 
SCS, where (E c — li) > 4Eg and quantum effects limit 
the role of short range potential. To model such a case we 
continue to use classical Hamiltonian Eq. ([3]), but trun- 
cate the Coulomb potential to V(r) — e 2 / K\r 2 + a^) 1 / 2 . 
Note that Eq. (|3| does not include the kinetic energy of 
electrons and holes in conduction and valence bands and, 
therefore, aims only at description of the low temperature 
(k B T < Eg) physics of SCS. 

Below we use dimensionless units, for r, ag, H, E g , 

— 1/3 

ksT measuring all distances in units of N D and all 
energies in units of e 2 Np 3 /n. Thus, Eq. |3| now can be 
understood as dimensionless, where E g 3> 1 and V(r) — 
(r 2 + o|) -1 /2. For TI with the energy gap 0.3 eV, K = 30 
and N D = 10 19 cm" 3 , N D 1/3 = 4.6 nm, e 2 N^ /3 /K = 10 



meV, so that dimensionless E g = 30. We could not model 
Eg = 30, because in this case the very large correlation 
length of long range potential, R, leads to large size effect. 
Instead, we run more modest E g — 15, for which the size 
effect requires extrapolation only for K = ffl Our goal is 
to find activation energy A and estmate Tes as functions 
of K or //. 

We search for the set {fii, fi} that minimizes H and 
use such a set to calculate the DOS and the conductivity. 
We start from the neutral system of all populated by elec- 
trons (negatively charged) acceptors (rij = l,<fc = — 1), 
of equal number of randomly chosen 20000-ftT empty 
(positively charged) donors (rij = 0, % = 1), and of 
20000(1 - K) filled (neutral) donors (rij = l.ft = 0). 
Charged donors and acceptors create a random poten- 
tial whose magnitude exceeds E g . In order to screen the 
Coulomb potential fluctuations, some electrons leave ac- 
ceptors for donors. At any stage of this process there are 
two types of occupied states - neutral donors and nega- 
tively charged acceptors, and two types of empty states - 
positively charged donors and neutral acceptors, respec- 
tively. Electrons may hop from an occupied impurity to 
an empty one. If the proposed move lowers the total sys- 
tem energy H, then it is accepted, otherwise it is rejected. 
To check whether H goes down, for a given electron oc- 
cupation numbers it is convenient to introduce 
the single-electron energy state, Ei, at a given impurity i 

For all i, j with rij = 1 and rij = we check that ES 
pseudo-ground state stability criterion is satisfied: 

e j -s i -V{r ij )>0. (5) 

If this criterion is not satisfied we move electron from im- 
purity i to j and recalculate all This process is done by 
looping all possible pairs of impurities i, j with rij = 1 and 
rij = and is continuted until no single-electron trans- 
fers are possible that lower H. The final arrangement of 
electrons can be called a pseudo-ground state, because 
higher stability criteria of ground state are not checked. 
Pseudo-ground states are kn own t o describe properties of 
real ground state pretty wel l 18 * 24 l Results below are ob- 
tained at E g = 15, a B = 1 for K = l, 0.99, 0.98, 0.97, 0.96 
and 0.95 and are averaged over 100 realizations of impu- 
rities coordinates. 

For a pseudo-ground state we find Fermi energy /iasa 
half distance between minimum empty and maximum oc- 
cupied energy e. Fig. [2] shows how the Fermi level ll{K) 
shifts from the middle of the gap towards the conduction 
band bottom with growing 1 — K . At 1 — K > 0.01 this 
dependence is in reasonable agreement with prediction 
of single band theory (the theory which ignores valence 
band and acceptors) ^ that E c - Li = A(l - Ky 1 / 3 . 
Note, however that for heavily doped SCS the coefficient 
Ah ~ 1.4 is twice smaller than the coefficient A[ ~ 2.8 ob- 
tained in Ref.^ for a lightly doped SCS, where qb <§C 1. 
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FIG. 2. (Color online) Fermi level /i as a function of 1 — K 
for a_g = 1 and i? 9 = 15. The size of dots characterizes the 
uncertainty. 
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FIG. 3. (Color online) Dimensionless single-electron DOS g(e) 
in units [(1 + K)N/(e 2 N 1 ^ 3 /k)] as a function of e calculated 
from the Fermi level for as = 1 and E g — 15 at K = 0.95 
(blue) and 1 (red). Impurity states with e < are occupied 
and with e > are empty. At K = 1 the total DOS of 
impurities has donor-acceptor symmetry, which is lost with 
growing 1 — K. 



In this case the short range Coulomb interaction at dis- 
tances r <C N^ 1 / 3 leads to an additional contribution to 
/x of the same order of magnitude. 

To confirm our understanding of results for 1 — K > 
0.01 we obtained same results for the position of Fermi 
level n (and DOS of donors and conductivity, see below ) 
running a simplified one band model where all acceptors 
are assumed to be negative. Such program is similar to 
the classical impurity band program used in Chapter 14 
of Ref. 18 , but uses the redefined V(r). 

The resulting DOS of impurities is shown in Fig. [3] for 
K = 1 and K = 0.95. g*(e) is the DOS in the units of 
(1 + K)N]j/(e 2 Np 3 / k) and is normalized to unity. 

At K = 1, almost constant symmetric DOS between 



FIG. 4. (Color online) Dimensionless DOS g(e) for neutral 
(occupied by electrons) donors at e < and neutral (empty) 
acceptors at e > as a function of energy e calculated from 



the Fermi level for as = 
(green) and 0.98 (blue). 



1 and Eg = 15 at K = 1 (red), 0.! 



tribution of random potential from —E g /2 to E g /2 and 
of band edges E c and E v between to E g and between 
to — Eg correspondingly (see Fig. [Ik ). In the middle 
(at the Fermi level) one sees the ES (Joulomb gap^I. 

At K < 1 the DOS of impurities loses donor- acceptor 
symmetry it has at K — 1. As mentioned in introduction 
(see Fig. [T]) with growing 1—K acceptors become all filled 
and disengaged from screening. Acceptor DOS (leftmost 
peak) splits from the donor one, which in turn has two 
peaks separated by the Fermi level. The large right peak 
belongs to empty donors, while the small and narrow left 
peak belongs to occupied donors. The donor peaks are 
separated by the ES Coulomb gap. 

Growing with 1 — K disengagement of acceptors from 
screening is also illustrated in Fig. [4j where we show DOS 
g(e) for neutral donors and acceptors. If at K = 1 the 
total number of electrons and holes in puddles are equal, 
with the growing 1 — K the total number of electrons in 
electron paddles grows, while the total number of holes 
in hole puddles decreases. Thus, at 1 — K > 0.02 valence 
band practically plays no role in screening. 



III. NUMERICAL MODELLING OF HOPPING 
CONDUCTIVITY 



Once the energies {ei} are known, we evaluate the re- 
sistivity using the approach of the Miller-Abrahams re- 
sistor network^. Each pair of impurities i,j is said to 
be connected by the resistance Rij = Rp exp[2r^/^ + 
^ij/ksT], where the activation energy Ey is defined^ as 
follows: 
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FIG. 5. (Color online) The temperature dependence of the 
resistivity in the whole temperature range 200 > T* > 
0.03. The dimensionless resistance (lnp)* is plotted against 
(T*)~ 1//2 to illustrate that the resistivity follows the ES law 
at low temperatures. The dashed lines are the best linear fits. 



FIG. 6. (Color online) The temperature dependence of the 
resistivity in the high temperature range 200 > T* > 1. The 
dimensionless resistance (In p)* is plotted against (T*) -1 to il- 
lustrate that the resistivity is activated at high temperatures. 
The dashed lines are the best linear fits. 



The resistivity of the system as a whole is found using a 
percolation approach. Specifically, we find the minimum 
value R c such that if all resistances Rij with Rij < R c 
are left intact, while others are eliminated (replaced with 
R — oo), then there exists a percolation pathway con- 
necting opposite faces of the simulation volume. The 
system resistivity p(T) is defined as RcN" 1 ^ 3 . Here we 
concentrate on exponential term of resistivity p ignoring 
details of the prefactoi^. 

For K=0.95, 0.97, 0.98 and 1 at a B = 1 and E g = 15, 
the computed dependence of (lnp)* = (£/2) In R c /Ro 
is shown function of (T*) 1 / 2 in the huge range 

of temperatures 200 > T* > 0.03 in Fig. [FJ Here 
T* = 2kgT/t; is yet another dimensionless temperature. 
These notations are introduced to exclude any explicit 
dependence on ^. It is seen that at low temperatures 
0.3 > T* > 0.03, resistivity is well described by ES law 
Eq. pi) (with C ~ 4.4 at K = 1). The higher temperature 
range 200 > T* > 1 is plotted separately as a function of 
1/T* in Fig. [6]o. We find two activated regimes of hop- 
ping conductivity. At high temperatures 200 > T* > 50 
we see the large activation energy E a ~ (E c — p) and 
at intermediate temperatures range E g > T* > 1 we see 
much smaller activation energy A = 0.3(£' c — p). 

The first activation energy E a does not have any 
physics meaning for a real SCS, because at k B T > E g 
conductance of SCS is actually not due to hopping, but 
due to high energy free carriers, which are not taken into 
account by energy Eq. Q (see RefP). In contrary to 
E a , the second activation energy A = 0.3(_E C — p) makes 
full physics sense and should be seen in an real experi- 
ment. The origin of this activation energy for the hopping 
transport is also explained in Chapter 8 of RefP^). At 
T <C Eg electrons optimize their conductivity using for 
hopping impurities energetically close to the Fermi level. 
Eventually at very low temperatures such opitmization 



leads to ES conductivity. However, when donor energies 
are slowly modulated by the long range potential there 
are large areas which do not have donors with energies 
close to the Fermi level and tunneling through them is 
slow. Therefore, there is a range of temperatures where, 
electrons use only nearest neighbor donors for hopping, 
while activating to donors located at the percolation level 
of nearest neighbor percolation. We find activation en- 
ergy from the Fermi level to the nearest neighbor perco- 
lation level studying hopping activation energy A. In a 
heavily doped semiconductor this energy is not descrim- 
inable from activation energy of electrons from the Fermi 
level to the conduction band percolation level E e . (Of 
course, holes are activated from the Fermi level to their 
percolation Eh as well so that A = 0.3(p — Eh)). 

We verified that hopping conduction modelling cor- 
rectly predicts activation energy of the band transport 
by direct of calculation of the percolation level E e . For 
this purpose we created a cubic lattice with a small lattice 
constant A r_1 / 3 /3. At every site of this lattice we calcu- 
lated potential of all charged impurities and then found 
lowest energy E e at which percolation over this lattice 
takes place. Activation energy of the band transport was 
again close to A = Q.3(E e — p). This result is also close to 
one obtained in Ref.^ based on estimate of percolation 
level for a generic random long range potential. 

In Fig. [7]we plot the function A(E c — p,) for all the val- 
ues p(K) obtained for K = 1,0.99,0.98,0.97,0.96,0.95 
we studied. We see that equation A ~ 0.3(£' c — p.) holds 
well for all K in this interval. 

So far we emphasized results, which do not explicitly 
depend on £. Actually, a magnitude of £ is necessary to 
calculate Tes- We argue now that in a TI £ is quite large 
leading to prominent role of VRH. If electron with energy 
close to the Fermi level were tunneling from an electron 
puddle to a distant one along the straight line it would 
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FIG. 7. (Color online) The activation energy A at if = 
1.0,0.99,0.98,0.97,0.96, and 0.95 (from right to left). The 
dashed line is the best linear fit A ~ 0.3(E C — fi) 



tunnel through high barriers and its wave function would 
decay with £ -C a B - Actually, a tunneling electron can 
use the same geometrical path as a classical percolating 
electron with energy A above the Fermi level which avoid 
large barriers. We assume that along such a path tun- 
neling barriers V are uniformly distributed in the range 
< V < A and neglect contribution of curvature of this 
path into action. Integration over V then gives (from here 
on we return to normal units) £ = %/ (StoA/9) 1 / 2 and 
k B T ES = 4.2(e 2 /Kfi)(mA) 1 / 2 . For a TI with a B = N D 1/3 
we get T ES = 4.2[(e 2 iV^ /3 /K)A] 1 / 2 . For A varying be- 
tween 1 and 2.5 e 2 N^ 3 /k as shown in Fig. [^T^s changes 
from 4.2 to 6.6 e^N^/n. For k = 30, N D = 10 19 cto- 3 

e 2 N^ /3 /Kk B ^100 K and T ES varies from 420 to 1050 
K. In order to study VRH in TI samples experimentally 
one has to deal with large enough samples, where surface 
conductance is smaller than the bulk one 26 . 



IV. THERMOPOWER 

In the recent papei^ the authors studied activation 
energy of the bulk resistivity of series of samples of 
Bi2Te3_ a; Sea ; with different x and, therefore, different po- 
sitions of the Fermi level in the TI gap. They found that 
when the Fermi level sinks into the gap the activation en- 
ergy of resistivity A grows and reaches maximum at 40 
meV and then decreases. The increase of the activation 
energy A on both sides of the maximum is accompanied 
by the increase of the absolute value of the thermopower 
S. However, near the maximum of A the thermopower 
abruptly changes its sign. These findings in agreement 
with what one can expect when a semiconductor goes 
through the point of complete compensation. Here we 
would like to concentrate on the maximum absolute value 



of the thermopower, for example, at n-type side of the 
maximum. 

It is known that for flat bands n-type semiconductor 
with the Fermi level /i inside its gap the thermopower 
S = A/eT , where the activation energy A = E c — /i. For 
bended bands of a strongly compensated rt-type semicon- 
ductor one could think that S = A/eT, where the acti- 
vation energy A = E e — fi is determined by activation 
to percolation level E e . Actually it was arguecP^MU that 
the Peltier energy (heat) n = eTS is determined by the 
average potential energy of electrons (conduction band 
bottom) along most conducting one-dimensional perco- 
lation paths, n =< E — fi >. (We call a percolation path 
any line, where potential energy of electron is smaller 
than E e and we call a set of least resitive of these paths, 
which carry most of the current, as most conducting per- 
colation paths.) Thermopower of an open circuit follow- 
ing an individual percolation path can be obtained by 
integrating E — fx along this path. Among two paral- 
lel paths connecting points A and B, the more resistive 
one has somewhat larger open circuit thermopower and, 
therefore, drives circular current back through the least 
resitive one. This current reduces thermopower of the 
resistive path so that voltage between A and B is deter- 
mined by the more conducting path. 

If probability distribution of potential energy E on 
most conducting paths is the same as for unconditional 
probability distribution of E, which we call DOS g(E) 
above, we can use g(E) to calculate n and S. For ex- 
ample, in the case of a constant g(E) for /i < E < E e 
we get E s —< E — fi >= A/2 = (E e - fi)/2. This con- 
clusion was confirmed by the numerical experiment!^ for 
the case of a constant g(E). 

In a stronly compensated semiconductor one can use 
real g(E) found above. For example, at K = 0.95 one 
can use Fig. [2] to find that A = E e — pi ~ 3. Then using 
DOS shown at Fig. [3] one can check that the average 
energy in the range between in the interval < E < 3 
is < E — fi >~ A/2 = 1.5. Thus, our simple approxi- 
mate prediction is that n A/2 = E g /12. It is valid for 
all K < 0.98 we studied. In the future work we plan to 
verify this prediction by a numerical calculation of the 
thermopower as a function of K for potential created by 
donors and acceptors in strongly compensated semicon- 
ductor, which we generated in the modelling above. 

For the data of the papei!^ our prediction means that 
at T = 50K the largest thermopower S = H/eT observed 
should be of the order 25mVy50K= 0.5 mV/K in a good 
agreement with observed value S — 0.4 mV/K. 
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